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Abstract 

A model of interacting random walkers is presented and shown to give rise to pat- 
terns consisting in periodic arrangements of fluctuating particle clusters. The model 
represents biological individuals that die or reproduce at rates depending on the 
number of neighbors within a given distance. We evaluate the importance of the 
discrete and fluctuating character of this particle model on the pattern forming 
process. To this end, a deterministic mean-field description, including a linear sta- 
bility and a weakly nonlinear analysis, is given and compared with the particle 
model. The deterministic approach is shown to reproduce some of the features of 
the discrete description, in particular, the existence of a finite-wavelength insta- 
bility. Stochasticity in the particle dynamics, however, has strong effects in other 
important aspects such as the parameter values at which pattern formation occurs, 
or the nature of the homogeneous phase. 

Key words: Pattern formation, Fluctuations, Interacting particle systems, 
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1 Introduction 



Pattern formation in the presence of fluctuations [1,2] is a subject attracting 
much interest since the beginnings of instability studies. The most common 
way to model the process is in terms of stochastic field equations in which the 
deterministic part undergoes a pattern forming instability and a noise term 
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modifies the dynamics. In situations where fluctuations are of internal origin, 
as for example arising from the discrete nature of the particles integrating the 
system, a more fundamental description is in terms of the microscopic particle 
dynamics (an analogy between a particular competing individual system with 
a reaction-diffusion system is presented in [3]). Models of interacting particle 
systems leading to periodic pattern formation are however not abundant in the 
literature. We introduced recently [4] one such a model, inspired in the inter- 
actions of biological individuals that compete for the resources around them 
and generalizing a model introduced in the context of plankton populations 
[5] . It consists in a set of random walkers, each one reproducing or dying with 
probabilities depending on the number of other individuals in a neighborhood. 
In our previous work [4] we presented a detailed numerical study of the model 
in two dimensions whose main result can be summarized in the existence of 
a stationary periodic pattern of clusters of particles for a specific range of 
parameters. A continuum Langevin equation is then derived under certain ap- 
proximations for the density of particles, with a rather complicated noise term. 
A linear stability analysis is then performed on the deterministic (noiseless) 
equation that in two dimensions predicts the transition to periodic patters by 
a finite- wavelength instability. However, some numerical results, like the value 
of the parameters at which patterns emerge, are not well reproduced in the 
deterministic density equation. This reflects that the role of fluctuations (the 
noise term in the Langevin equation) maybe rather important depending on 
the properties of the system one wants to focus on. 

In this context, we present in this Paper additional results for the model, with 
more emphasis on the one- dimensional case, and compare it more carefully 
with descriptions in which fluctuations are absent: a mean-field approximation 
and its amplitude equation. The onedimensional case is relevant by itself (see 
for example a relevant experimental situations in [6] and references therein), 
and also because it allows for a simpler analytical treatment. The amplitude 
equation in this work is calculated for both the one- and two-dimensional 
situations, but it is only the onedimensional expression which turns out to be 
suitable for full comparisons with the numerics. It is also worth mentioning 
some other extensions with respect to the results shown in [4]: the analytical 
results are stated for a general Kernel (see Sect. IV) of the mean-field density 



equation, and two order parameters are identified in the model, v//(0) and 



yl{km), which inform, respectively, about the empty-active transition and 
the pattern formation transition (see Sect. III). 

The Paper is organized as follows: In Sect. 2 we present the particle model, 
and show numerical results on its behavior in Sect. 3. Section 4 deals with the 
mean-field approximation, including its linear and weakly nonlinear analysis. 
We close with the Conclusions. 
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2 A model of nonlocally interacting individuals 



We consider here both the one-dimensional situation {d — 1), in which the 
system, consisting of N{t) particles, is contained in a line, and the two- 
dimensional one {d = 2), in which the N{t) particles are in a square. In both 
cases wc use periodic boundary conditions. At the starting time t = an ini- 
tial population of N{0) particles is randomly distributed. At time t, when the 
population is N{t), one particle is selected at random (say particle j) and it 
dies (thus disappearing from the system) with probability pj, reproduces (i.e. 
it replicates itself) with probability qj, or is left unchanged with probability 
1 —pj — Qj- In the case of reproduction the newborn is placed at the same loca- 
tion as the parent particle. The process is repeated N{t) times. After this, each 
particle moves independently in a random direction for a distance drawn from 
a Gaussian distribution of standard deviation a, then time is incremented an 
amount r {t ^ t + r), and the algorithm repeats with the resulting N{t + r) 
individuals. The random motion leads to diffusion with a diffusion coefficient 
D = a"^ /2t. The key point is the election of the probabilities pj and qj or, 
equivalently, the death and reproduction rates (5j = pj/r and \j = qj/r, re- 
spectively. We choose a constant death rate Pj — /3o, Vj, but the reproduction 
rate of particle j depends of the number of particles A^^ within a distance it! 
from its position: 



Ao - gNi if K < Xo/g 
if A^^ > Xo/g 

This kind of interaction models a slowing down of birth rates (until total 

suppression) in regions of high particle density. This is rather appropriate to 
model biological populations that compete for food or other resources present 
in a neighborhood of their positions, or that release toxic chemicals. 

We choose Xq + (Sq = 1/r, and introduce the maximum net growth rate /i = 
^0 ~ Po- For fixed r, the value of /i fixes both Aq and Pq. Since pj and qj 
are probabilities, we see that /xr G [—1, 1]. One can estimate an equilibrium 
average density in this model by imposing that at each site birth and death 
compensate on average: Aq — g{Nj^) — Pq. By assuming a uniform average 
density, 0, one can express the number of neighbors as (Nj/) = (phR^ {h is 2 
in one dimension and vr in two dimensions, so that hR^ is the length or the 
area of the neighborhood of a particle; we have not discounted the central 
particle from the number Nj^, since a difference of one particle is irrelevant in 
the regime this argument is intended to describe). Thus 




But this estimation completely neglects inhomogeneities and diffusion. Fluc- 
tuations are known to produce inhomogeneities that alter estimations such as 
(2). This is specially true in situations like this where, in addition to (2), there 
is another homogeneous state where birth and death processes also compen- 
sate: the empty state = {N{t) — 0), for which both rates vanish. This 
is an absorbing state [7] since further evolution is impossible if the system 
reaches it. A phase transition to the absorbing state is expected by reducing 
the net growth rate /j, from a higher value. The character of the transition and 
the location of the transition point are not well reproduced if fluctuations are 
neglected [7,8]. 



3 Numerical observations 

We first consider the one-dimensional case. The particle model described above 
is run on a ring of length L = 1. For smaller than a critical value fiQ, ah initial 
particle distributions that we have checked end up attracted by the absorbing 
empty state. An illustrative way of showing one- dimensional particle dynamics 
consists in plotting the particle positions in the horizontal axis and the time 
evolution of these in the vertical. Proceeding in this way. Fig. 1 shows an 
example of the dynamics of the emptying process, strongly reminiscent of 
directed percolation [8] below threshold. 



By increasing the net growth, fi, two situations appear: If the diffusion coeffi- 
cient is large enough, an active phase, roughly homogeneous as shown in Fig. 
2, appears above a critical value hq. The value of iio is larger than the one 
required to equilibrate maximum birth rate and deaths {/i — 0). Moreover the 
average density is smaller than the predicted by the crude estimation (2). 



The other situation appears for smaller diffusion coefficients: now increasing 
/X above fiQ > particles organize in clusters which are roughly equidistant, so 
that a periodic pattern results (Figs. 3 and 4). The configurations are rather 
noisy, specially close to the transition (Fig. 3), but a periodicity in the pattern 
is clear. 



In the two-dimensional case the general scenario is similar, and was already 
described in [4]. Patterns consist in fluctuating particle clusters arranged in a 
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X 



Fig. 1. Time evolution (in vertical) of the particle positions (in horizontal) for 
fj, = 0.4, D = 10^^, R = 0.1, and g = 1/50, starting from an initial population 
of 1000 particles. Units of time are such that r = 1, and system size is L = 1. The 
final state is empty of particles 




X 



Fig. 2. Idem as in Fig. (1) except for n = 0.46 and D = 3 • 10 ^. An active homo- 
geneous phase has developed. 
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Fig. 3. Idem as in Fig. (1) except for /x = 0.50. Noisy pattern formation has occurred. 




Fig. 4. Idem as in Fig. (1) except for jj, = 0.90. A rather well defined pattern is 
reached soon. 

hexagonal lattice. We show a snapshot in Figure 5. 
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Fig. 5. A snapshot of the fluctuating two-dimensional pattern of particle clusters, 
obtained in a square of side L = 1 for = 0.7, D = 10~^, i? = 0.1, 5 = 1/50, and 
r = 1. 

In order to be more quantitative, we introduce a structure function 



m 



N{t) 



(3) 



The sum is over particles, at positions Xj, and the average is a temporal average 
on a long-time state, performed to improve statistics. Maxima in I{k) identify 
the wavcnumbcrs, k = |k|, associated to periodicities in the system. In any 
nonempty state there is also a peak in k = giving the square of the mean 
particle density. The values 7(0) and I{kM), where /cm is the nonvanishing 
wavenumber at which there is a maximum, act as two order parameters: the 
first one for the transition from the empty to the active phase, and the second 
for the pattern formation. For Figs. 3 and 4, kM = 147r, corresponding to 
the observation that the pattern consists of 7 cells (a pattern of n cells is 
characterized by a wavenumber kM — 2nn/L). In Fig. 6 we plot the value 
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Fig. 6. Values of \/I{0) (squares) and \/I{kM) (circles) as a function of fx for the 
one-dimensional particle model aX D = 9.2 • 10^'', R = 0.13, g = 1/120, r = 1 and 
L = 1, for which = IOtt. The upper thick line is (j)s, from (32), and the lower 
one is 0^, from (33). 

of \J and ^ I{kM) as a function of /i in a different case for which ku — 
IOtt (i.e. 5 cells). Both order parameters become nonvanishing at the same 
absorbing transition, although with rather different amplitudes and behavior. 



4 Mean-field description 



We try now to understand the pattern forming processes described in the pre- 
vious Section. It is obvious that the patterns are rather noisy, suggesting that 
fluctuations play a role in its development. It is known from related models 
[9,5,10] that, in addition to its impact on the absorbing transition, reproduc- 
tive fluctuations have another important effect arising from the asymmetry 
between birth and death: Birth is a multiplicative process that increments the 
density in regions where it is already high, whereas death can occur anywhere. 
Thus, even if the average rates are equal, under reproductive fluctuations there 
is a tendency of the particles to organize into clusters of high density, lead- 
ing to very inhomogeneous configurations when diffusion is not strong and is 
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the only homogenizing force [9.5.10]. This inhomogeneity is amplified when 
birth rate exceeds death, but it is not observed in the common situations in 
which the growth is saturated by local interactions [8] . Thus, it is not obvious 
whether the pattern development in our model arises from the stochastic mi- 
croscopic fluctuations, or if rather it is caused by some deterministic instability 
arising from the nonlocal interactions. 

To address this question we can write down a mean-field-like description of 
the model, which completely neglects fluctuations, and check if the instability 
is present there. The mean-field equation is written in terms of an expected 
density 0(x, t) as follows 

at</.(x, t) = A^0(x, t) + DV'cj>{^, t) - gcj>{^, t) I ^(x - y)0(y, t)dy . (4) 

s 



The first and second terms are the standard mean-field description of the net 
growth and the diffusion processes. The last one is a nonlocal contribution 
associated to the saturation produced by the neighborhood interactions in 
our model if the kernel G(x) is given by 



G(x) = 



1 if Ixl < R 

(5) 

if X > R. 



S is the domain were the system is defined. Equation (4) was more system- 
atically deduced in [4] from our microscopic model, as the deterministic part 
of a Langevin equation containing rather complex noise terms. Equation (4) 
and related ones have also been proposed independently to model directly the 
macroscopic behavior of a variety of biological systems [11,12,13,14]. In these 
cases, other kernels were considered in addition to (5), and for some of them 
pattern formation was observed. In the following we establish properties of (4) 
for general interaction kernels such that G = G(|x|), although the final results 
will be only commented for the step kernel given by (5). We first estabhsh 
general relationships, then perform a linear stability analysis of a homoge- 
neous solution, and then consider the weakly nonlinear behavior. We restrict 
to domains S either unbounded or with periodic boundary conditions. 



4-1 General relationships and homogeneous solutions 



It is useful to write Eq. (4) in non-dimensional form to have a clearer idea of 
the independent parameters involved. We rescale the units of space, time, and 
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(f) as follows: 



In this way: 

at*(u, s) = /i*(u, t) + V^^(u, t) - ^(u, t) J H{n - v)^(u, t)dy . (7) 

E 

E is the rescaled domain, and H{vl) is the rescaled kernel: H{\i) = G(x = Ru). 
In this mean-field description, The only relevant parameter turns out to be 

For future reference we introduce 

H{q) = J H{u)e''^ ''du , (9) 

E 



and the notation h = H{0) = H{u)du. Because of the symmetry assumed 
for G'(x), H{q) depends only on the modulus of the wavenumber g = |q|. For 
our step kernel we have, in the one-dimensional case: 

^(g) = 2^, h^2 (10) 



and for the two-dimensional case 

H{q) = 27r^^ , h ^ tt . (11) 



Ji(g) is the first order Bessel function. 

The trivial solution ^^=0 of (7) becomes unstable for jl > 0. Then a uniform 
solution appears: 

*(u,i) = ^ = ^, (12) 



corresponding to 
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This reproduces the ciTide approximation (2), performed for the particle model. 
It is convenient to introduce the deviation -0 with respect to the homogeneous 
solution: 



^(u,t) = ^+^(u,t), (14) 



so that (7) reads 

dMn,t)^Vl^{u,t) 

-u J H{u - v)^/'(v, t)dv - iIj{u, t) J H{u - v)^(v, t)civ . (15) 

E E 

4-2 Linear stability analysis 



The linear stability analysis of the homogeneous solution can now readily 
performed by introducing perturbations of the form ip ~ exp(iq • u + X{k)t) 
and neglecting nonlinear terms. The perturbation growth rate is 

X{q)^-q'-uH{q) . (16) 



It is seen that a necessary condition for instability is that H{q) takes negative 
values (note that Eq. (16) generalizes to an arbitrary kernel the expression 
shown in [4] for the step kernel) . It is also sufficient, since instability will 
always be reached in this case by increasing sufficiently u (or jl). Thus, if 
H{q) < for some range of q, there is a critical value of (corresponding to 
a critical such that an instability to pattern formation occurs [12,4]. Since 
both (10) and (11) take negative values, we see that the mean-field model with 
the step kernel undergoes a pattern forming instability without the need of any 
stochastic noise, as already reported [11,4,12]. Figure 7 (left) plots expression 
(16) in the one-dimensional case, showing the change of sign of A(g) around a 
critical qc by increasing \x — vh. Figure 8 (left) shows a long-time numerical 
solution of (4)- (5) in one dimension, displaying a well-developed steady pat- 
tern. The right part of Fig. 8 shows the analogous structure in two dimensions: 
the instability leads to the formation of a pattern of hexagonal symmetry, as 
in the particle model. The corresponding two-dimensional dispersion relation 
(16) is shown in the right part of Fig. 7. 
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Fig. 7. The linear growth rate (16) of perturbations of wavenumber q to the homo- 
geneous solution (12). Left: d = 1. Right: d = 2. In both cases positive growth, i.e. 
instability, arises by increasing fl. 
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Fig. 8. Steady periodic patterns obtained by long-time integration of the mean-field 
model (4) in d = 1 (left) and d = 2 (right). In both cases fi = 0.7, D = 10-^ R = l, 
and g = 1/50. In the second case higher values of (j) are represented by lighter gray 
levels. 
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Eq. (16) explains also the absence of instability for other kinds of kernels, such 
as the Gaussian one, for which H{q) is positive definite [13,11,12] (we stress 
however that patterns may appear under front propagation, or in response 
to spatial inhomogeneities [13]). To make a more quantitative analysis, the 
value of Vc and the fastest growing wavenumber Qc are obtained from the 
simultaneous solution of 

2q, + u,H'(q,)^0. (17) 

Numerical solution of these equations leads, for the one-dimensional case and 
step kernel: 

gc ~ 4.078 , z/c ~ 42.10 , (18) 

so that jjLc ~ 84.20. In two dimensions: 

qc ^ 4.779 , i/c ~ 58.948 , (19) 



and jlc ~ 185.192. Prom (6), the unsealed values k of the wavenumbers are 
given by k = q/i?. 

These predictions can be compared with the numerical observations of the 
particle model. First, the maximum of the structure function for the patterns 
in the Figs. 3 and 4 was &t k — lAn 43.98. This is precisely the first 
wavenumber compatible with the periodic boundary conditions immediately 
above the fastest growing mode in the mean-field equation for R = 0.1: kc = 
qc/R ~ 40.78. For the configurations of periodic character among the ones 
used to construct Fig. 6, the maximum in the structure function was at A; = 
IOtt 31.4159, to be compared with kc = qc/R ^ 31.369 for R = 0.13. 
The agreement is thus excellent in both cases. We can conclude that the 
essential mechanism of the pattern forming instability is well captured by the 
deterministic description, being the period of the pattern even quantitatively 
reproduced. For the two-dimensional case, the same agreement was already 
observed [4]. 

The situation is different for the location of the point of pattern onset. As 
for the location of the absorbing transition, the mean-field prediction and the 
particle model disagree. Typically, the transition from the empty state to the 
pattern is observed directly, without the intermediate of the homogeneous 
distribution. Only if diffusion is greatly increased the homogeneous state is 
observed (Fig. 2) but then the transition to patterns predicted by (18) occurs 
at 11 exceeding the maximum of the physical range /it G [—1, 1] (and it is 
indeed not observed in the particle model, which is a kind of agreement). 
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The reason for this disagreement, despite containing the correct instabihty 
mechanism, is easy to understand: In cases such as the ones in Sect. 3 with 
smaU D {D = 10~^) it happens that the predicted value of /ic (= 0.084) is 
much smaller than the //q of the absorbing transition (which is 0.45). Thus 
it is not strange than, as soon as the active phase becomes preferred above 
fiQ, it is aheady of the periodic type. The failures of the mean-field approach 
to fully describe the pattern behavior are then due to the strong impact of 
microscopic fluctuations on the absorbing transition. A second case is the one 
corresponding to Fig. 6. Here fic — 0.458 > /^o ~ 0.3. But what is observed 
is that, even when the system is below threshold, microscopic fluctuations 
excite the less-damped modes [1,2], so that a small amplitude noisy pattern is 
always present above the absorbing transition, as seen in Fig. 6. Fluctuation 
strength is relatively large here because of the low amplitude of the active state 
close to /ic, caused by the proximity to /Iq. Pattern formation is not a sharp 
bifurcation in this case, but a continuous crossover. Finally, in the situation 
in which the absorbing transition and the pattern forming instability would 
be well separated, this last one falls outside the physical range /it G [—1, 1]. 

4-3 Weakly nonlinear amplitude equations 

The advantage of having a continuum model as (4) is the possibility to use 
the tools of pattern formation theory [15,16] to obtain analytic results. A 

first example has been the linear stability analysis performed before. We now 
estimate the pattern amplitude close to threshold by calculating the pertinent 
amplitude equations. 

To this end we start from Eq. (15), and perform the standard expansions in a 
small parameter e related to the distance to threshold: 



V'(x, t) = e^i(u, Ti, T2, ...) + e2V'2(u, t, Ti, T^, ...) + e^iPs{xu, t, Ti, T^, ...) + .. 



We have introduced two slow time scales: Ti = et, T2 = e^t. For simplicity we 
do not introduce additional space variables, needed when considering large- 
scale pattern modulations. 

To first order in e, one finds the equation 



Ti is the integral operator acting as 7if{\x) = J^H{u — v)/(v)(iv. Eq. (22) 
implies that ipi is a combination of periodic patterns of the critical periodicity 



U ^ + 61/1 + e^V2 + ... 



(21) 



L,^i(u, Ti, T2) = , with L, = dt-Vl + v^n . 



(22) 
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27r/?c: 

V^i(u,Ti,T2) = ^ A,(Ti,T2)e^^^ " , (23) 

3 

with IQjI = qc- 

At orders and we find 

= -<9ri^i - yiU^i - V'l'Wi/'i (24) 

= -9ti'02 - ^TaV'l " " ^'2'W'01 " - ■0l'^'02 ■ (25) 

We particularize to the one-dimensional case, in which the only possible plan- 
form of the form (23) is 

V'i(m,Ti,T2) = /l(Ti,T2)e'9=" + ^*(Ti,T2)e-*^^" . (26) 

A* is the complex conjugate of A. Substituting in (24) one finds, as usual in 
one-dimensional patterns, that boundedness of the solutions implies i/i — 0, so 
that e oc (i/ - UcY^^ and A{Ti,T2) = A{T2). By solving (24) and substituting 
in (25), one finds the condition to avoid resonant terms: 

Ot.A = -U2H(qc)A - k\A\'^A , (27) 

which is the desired amplitude equation. Here, 

K = H{q,) {Q + S) + H{2q,)Q + H{0)S 

-Hiqc) ^ -2H(qc) , , 

^ql^v,E{2q,) u,H{'^) 

All these quantities are pure non-dimensional numbers. Their values, in ad- 
dition to the ones given in (18), are H{Qi) — h — 2, H{qc) ~ —0.3950, 
H{2q^) ^ 0.2341, Q = 5.172 • lO'^, S = 9.383 • 10-^ and k = 0.01423. 
The content of (27) is more clearly seen by reintroducing the original time 
variable t — e~^T2 and the amplitude B{t) — eA{T2) so that 

B{t)^{u-u,)[-H{q,))B-K\B\''B (29) 

The equation predicts that instability develops for v > Uc and leads to a 
pattern of arbitrary phase, and amplitude 



|5| = V("--c)^^ (30) 
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that has bifurcated supercritically from the homogeneous sohition (12). Re- 
turning back to dimensional variables via (6), the periodic steady solution for 
the original density can be reconstructed as 

<l>{x, = 0, + 0p (e^(^=-+'^) + e-^(*^=-+'^)) + - ^^,) . (31) 

with 



2gR 
and 



(32) 



0p = v/^^V^- (33) 



99 is an arbitrary phase. We can compare this analytic prediction with the 
particle patterns. The structure function (3) can be interpreted as the power 
spectrum of a density of the form Ya=\ (^(x — x^). Comparing with the power 
spectrum of the mean- field density (31), we see that should be related 

to 0s, and ^jliku) to (pp. We remind that /cm was well reproduced by the 
mean-field value of kc- Comparison is performed in Fig. 6. We see that the 
agreement in the amplitudes is poor (even if it was expected to be valid only 
close to lie) although some of the trends are reproduced. The disagreement is 
linked to the failure to capture the correct transition points. 

The amplitude equation can also be calculated in 0? = 2, by using an hexagonal 
planform [15,16] instead of (26): 

V'i(u,ri,r2)= Yl MT,,T2)e'^-^ + C.C. (34) 

J=l,2,3 



{j = 1, 2, 3), are three wavenumbers of modulus Qc at 120 degrees from each 
other, and c.c. means complex conjugate. The problem with such calculation 
is that generically the bifurcation to hexagons is subcritical, so that there is no 
guarantee that the amplitude equation will capture the saturation that occurs 
far from the instability point. By using (24)-(25) one finds, 

Bi{t) = {u- Ue) {-H{qc)) B, + 2 {-H{qc)) B^Bs + 0{\Bj\'Bi) , (35) 

and two other equations obtained by permuting cyclicly the indices 1, 2, 3. As 
in (29), Bi is the coefficient of the Fourier mode of wavenumber q, in the 
expansion of t/'(x, t). Eq. (11) should be used for H{qc)- The sign of the third 
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order terms not written explicitly in (35) is such that they do not saturate the 
pattern. Thus, expansion should be followed at higher orders and (35) is of not 
much utility. Nevertheless, the presence of the second order term confirms that 
the bifurcation in the mean-field model is subcritical, and the positive sign of 
its coefficient (remember that H{qc) < 0) implies that the bifurcating solution 
leads to HO or positive hexagons [16] (i.e. points of high density forming an 
hexagonal lattice), as observed in the particle model. 



5 Conclusion 



In this work we have presented new results on a model recently introduced 
by the authors [4] which considers particles reproducing and dying at rates 
depending on the number of individuals that every particle has in its i?-ranged 
neighborhood. The study has been performed at the level of the density equa- 
tion derived for this model, and mainly for the one-dimensional case, since 
analytical results are simpler here. 

As the general conclusion of the comparison between the particle and the 
mean-field model we can say that, although the instability appears to be 
deterministic in origin and the mean-field approach captures correctly the 
wavenumber, fluctuations are important and affect transition points and am- 
plitudes. This limits the usefulness of standard pattern forming theory tools, 
usually developed for continuum models as in the case of weakly nonlinear 
analysis, and ask for further studies in which noise terms are explicitly consid- 
ered [4]. Fluctuations may have also another interesting effect: even in cases 
in which there is no deterministic instability, as for example when a Gaus- 
sian kernel is used in (4) [13,11,12], it is very likely that a corresponding 
particle model would display a small amplitude noisy pattern, arising by the 
less-damped models excited by fluctuations. 
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